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1.  Introduction 


The  Thalmann  Exponential-Linear  Decompression  Algorithm  (EL-DCA)  was  developed 
with  the  VVal-18  and  HVal-21  parameter  sets  at  Navy  Experimental  Diving  Unit  (NEDU) 
to  produce  0.7  atmospheres  absolute  (ata)  constant  P02-in-N2  and  0.7  ata  constant  P02- 
in-He  decompression  tables.  Thalmann  published  FORTRAN  source  code  for  the  initial 
operating  version  of  the  EL-DCA  to  compute  decompression  tables  in  19831,2  and 
described  the  underlying  theory  in  1984. 3  Additional  FORTRAN  code  for  the  programs 
and  subroutines  of  the  EL-DCA  —  code  required  to  include  surfacing  repetitive  groups  in 
computed  decompression  tables  and  to  compute  surface/shallow  interval  credit  tables  to 
support  repetitive  and  multilevel  diving  —  was  described  and  published  still  later.4  5 

Thalmann’s  original  EL-DCA  was  extensively  modified  at  NEDU  to  develop  MK  16 
MOD  1  N2-02  and  He-02  decompression  tables6,7,8  and  Submarine  Rescue  and  Diving 
Recompression  System  Operator/Tender  decompression  tables.9  The  EL-DCA  was 
further  modified  and  used  with  a  modified  VVal-18  parameter  set,  designated  VVal-18M, 
to  produce  an  integrated  set  of  air  diving  decompression  tables  for  Revision  6  of  the 
U.S.  Navy  Diving  Manual. 10,11  This  report  describes  the  design  of  the  EL-DCA  software 
used  to  produce  the  integrated  air  decompression  tables  and  documents  the  software 
modifications  made  since  publication  of  Thalmann’s  original  reports. 

Principal  functional  enhancements  to  Thalmann’s  original  EL-DCA  include  support  of 
the  following: 

•  input  of  program  setup  parameters  from  a  user-editable  program  initialization  file, 

•  last  allowed  decompression  stop  depths  at  integer  multiples  of  the  stop  depth 
increment  other  than  one, 

•  in-water  oxygen  decompression  with  user-setable  inspired  oxygen  fraction  and 
automatic  insertion  of  air-breathing  breaks  (automatic  air/02  cycling), 

•  inclusion  of  travel  time  between  decompression  stops  in  tabulated  stop  times, 

•  rule-based  inclusion/exclusion  of  surfacing  repetitive  group  designators  with 
tabulated  schedules, 

•  automatic  depth-dependent  ascent  rates, 

•  calculation  of  schedules  with  instantaneous  ascent  to  surface  from  depths  other 
than  the  last  allowed  decompression  stop  or  without  instantaneous  ascent  to 
surface  from  any  depth, 

•  use  of  compartmental  saturation/desaturation  half-time  ratios  other  than  unity 
only  when  the  inspired  oxygen  fraction  is  greater  than  or  equal  to  a  threshold 
value  specified  in  the  program  initialization  file, 

•  support  for  threshold  inert  gas  tension  overpressure,  PBOVP,  different  at  surface 
from  that  at  depth, 
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•  computation  of  surface  interval  credit  tables  with  surfacing  letter  final  repetitive 
group  designators, 

•  computation  of  residual  nitrogen  times  for  surface  interval  credit/residual  nitrogen 
time  tables,  and 

•  output  of  all  tabulated  schedules  in  Augmented  NMRI  Standard  (ANS)  format 12 
for  processing  by  probabilistic  models  of  DCS  incidence  and  time  of  occurrence. 


The  Thalmann  Algorithm  decompression  table  generation  software  consists  of  two  main 
programs:  a  decompression  table  generator,  TBLP7R,  and  a  surface  interval  credit  and 
residual  nitrogen  time  table  generator,  SICT7.  A  third  program,  BLDCMPTBL,  was 
developed  to  supplement  TBLP7R  for  production  of  the  Air  Decompression  Table  in 
Revision  6  of  the  U.S.  Navy  Diving  Manual.  Each  of  these  programs  is  described  in 
separate  sections  of  this  document.  The  essential  features  of  program  design  and 
operation  are  depicted  in  information  and  process  flowcharts  with  only  minimal  narrative 
elaboration.  These  features  are  discerned  by  stepping  through  the  flowcharts  beginning 
at  the  starting  termination  of  the  TBLP7R  MAIN  chart  (Figure  1),  the  BLDCMPTBL 
MAIN  chart  (Figure  42),  or  the  SICT7  MAIN  chart  (Figure  47).  Supporting  procedures, 
for  which  separate  flowcharts  are  given  in  most  cases,  are  each  indicated  by  a 
procedure  name  in  a  double-sided  box: 


Procedure  Name 


The  flowcharts  for  these 


procedures  are  readily  found  by  reference  to  the  list  of  figures  on  pages  ii  and  Hi  of  this 
document.  An  annotated  list  of  program  variable  names  is  given  in  Appendix  A. 


Many  procedures  described  in  this  software  design  document  also  support  operation  of 
the  U.S.  Navy  Thalmann  Algorithm  Dive  Planner  (NDP),  sometimes  with  minor 
modifications.  Such  modifications  are  included  in  the  present  documentation  with  an 
“NDP  only”  designation  as  appropriate. 


In  both  TBLP7R  and  SICT7,  all  internal  calculations  are  completed  with  pressures  in 
units  of  feet  of  seawater  (fsw).  Surface  pressure  is  33  fsw  (SURFP_FSW  =  33)  in 
accord  with  the  value  used  in  Thalmann’s  original  code.  Internal  conversions  of  meters 
of  seawater  (msw)  to  fsw  are  made  with  a  conversion  factor  of  1 .0/0.3048  =  3.2808 
fsw/msw  (MSWpFSW  =  0.3048),  also  as  used  in  Thalmann’s  original  code.3  Maximum 
permissible  tissue  tensions  (MPTTs)  in  units  of  fsw  for  depths  in  msw  in  Thalmann’s 
original  VVal18.dat  file  were  also  computed  with  the  3.2808  fsw/msw  factor. 


2.  Decompression  Table  Generator  (TBLP7R) 

FORTRAN  code  for  the  original  version  of  TBLP7R  was  published  in  NEDU  TR  1-831 
with  addition  of  an  erratum  shortly  thereafter.2  Subroutines  to  augment  this  code  and 


a  Thalmann’s  3.2808  fsw/msw  conversion  factor  value  is  based  on  a  seawater  density  of  1 .02  gm/ml, 
apparently  rounded  from  the  value  of  1 .02480  gm/ml  given  by  the  Undersea  and  Hyperbaric  Medical 
Society  for  seawater  at  4°C.12  With  this  latter  unrounded  value,  the  conversion  factor  value  is  1 .0/0.30632 
=  3.26458  fsw/msw.  In  comparison,  surface  pressure  of  33  fsw  is  in  accord  with  the  value  of  33.0060  fsw 
obtained  with  the  seawater  density  of  1 .02704  gm/ml  reportedly  used  by  U.S.  submarines.14 
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support  computation  of  surfacing  repetitive  groups  and  residual  nitrogen  times  for 
repetitive  diving  were  described  in  NEDU  TR  9-85, 4  with  accompanying  code  in  NEDU 
TM  85-12. 5 

TBLP7R  processes  an  ASCII  input  file  of  instructions  for  the  schedules  that  are  to  be 
included  in  the  product  decompression  table.  This  file  must  be  in  Planner  Coordinate 
Format  (PCF)  described  in  Appendix  B.  Depending  on  the  content  of  the  PCF  input  file, 
TBLP7R  produces  either  a  table  of  decompression  schedules  or  a  Sub-NoD  table  of  no¬ 
stop  limits  with  corresponding  no-stop  times  required  to  reach  the  top  of  each  surfacing 
repetitive  group  at  each  dive  depth.  To  produce  a  table  of  decompression  schedules, 
the  PCF  input  file  should  contain  templates  for  the  schedules  to  be  computed  arranged 
in  order  of  increasing  bottom  time  in  groups  of  increasing  dive  depth.  To  compute  a 
Sub-NoD  table,  the  PCF  input  file  should  contain  a  template  for  one  dive  to  each  depth 
that  is  to  be  included  in  the  table.  Each  of  these  templates  should  then  prescribe  the 
compression  to  bottom  depth  with  a  dummy  bottom  time  at  that  depth  and  the  “SD” 
option  after  specification  of  the  gas  to  be  breathed  at  that  depth. 

TBLP7R  operation  starts  with  execution  of  the  procedures  illustrated  in  Figure  1 . 
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Figure  1.  TBLP7R  main. 

Program  switch  settings  and  other  parameters  input  by  the  Readjni  routine  from  a 
DCA_Params.ini  file  during  TBLP7R  startup  are  given  in  Table  1  with  brief  descriptions 
from  Gerth  and  Doolette.11 
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Table  1.  Parameters  input  by  Readjni  routine  during  TBLP7R  startup 


Parameter 

Description 

PAC02 

Arterial  C02  partial  pressure  (fsw) 

PH20 

Tissue  water  vapor  pressure  (fsw) 

PVC02 

Venous  C02  partial  pressure  (fsw) 

PV02 

Venous  02  partial  pressure  (fsw) 

AMBA02 

Ambient-arterial  02  partial  pressure  gradient  (fsw) 

PBOVP 

Threshold  compartmental  inert  gas  tension  overpressure 
(fsw)  with  respect  to  saturation  at  the  prevailing  ambient 
hydrostatic  pressure  (fsw)  for  transition  from  exponential 
to  linear  gas  exchange  kinetics 

sPBOVP 

PBOVP  at  surface  (fsw) 

LST_DOMode 

Last  stop  dropout  mode 

•  LST_DOMode  =  0:  In  last  decompression  stop  time 
calculations,  commands  allowance  for 
instantaneous  ascent  to  surface  from  a  depth  equal 
to  the  decompression  stop  depth  increment 

•  LST_DOMode  =  1 :  In  last  decompression  stop  time 
calculations,  commands  allowance  for 
instantaneous  ascent  to  surface  from  the  last 
allowed  decompression  stop  depth.  If  the  last 
allowed  stop  is  at  20  fsw  with  10  fsw  increments 
between  stops,  resultant  last  stop  times  are  longer 
than  those  obtained  with  LST_DOMode  =  0  by 
amounts  that  increase  with  decreasing  ascent  rate. 

•  LST_DOMode  =  2:  In  last  decompression  stop  time 
calculations,  commands  accounting  for  gas 
exchange  during  the  entire  ascent  from  the  last 
allowed  decompression  stop  depth  to  surface. 
Resultant  last  stop  times  are  shorter  than  those 
obtained  with  LST_DOMode  =  0  or  LST_DOMode  = 

1  by  amounts  that  increase  with  decreasing  ascent 
rate. 

SRF_CNTRLT_MODE 

Controlling  compartment  for  surfacing  mode 

•  SRF_CNTRLT_MODE  =  0:  The  controlling 
compartment  for  surfacing  is  the  highest  half  time 
(slowest)  compartment  with  gas  tension  in  excess  of 
its  surfacing  M-Val  when  the  diver  arrives  at  the  last 
decompression  stop  depth  during  ascent. 

•  SRF_CNTRLT_MODE  =  1 :  The  controlling 
compartment  for  surfacing  is  that  with  gas  tension 
closest  to  its  surfacing  M-Val  when  the  diver  leaves 
the  last  decompression  stop  depth  during  ascent. 

5 


Parameter 

Description 

RGD_SPRSS 

Repetitive  Group  Designator  (RGD)  suppression  mode 

•  RGD_SPPRSS  =  0:  Suppress  RGDs  within  a  dive 
depth  group  for  all  schedules  with  bottom  times 
longer  than  the  first  with  a  reference  compartment 
inert  gas  tension  on  leaving  bottom  greater  than  or 
equal  to  60  fsw  (bottom  number  repetitive  group  > 

60). 

•  RGD_SPPRSS  =  1 :  Suppress  RGDs  within  a  given 
dive  depth  group  for  all  schedules  with  bottom  times 
equal  to  or  longer  than  the  first  with  a  computed 

RGD  less  than  that  of  the  preceding  schedule. 

•  RGD_SPPRSS  =  2:  Suppress  RGDs  when  the 
controlling  compartment  for  surfacing  (see 
SRF_CNTRLT_MODE)  has  a  halftime  greater  than 
that  of  the  reference  compartment. 

BTMAX 

Maximum  Bottom  Time  (min) 

TATMAX 

Maximum  Total  Ascent  Time  (min) 

STIME 

Default  gas  switch  time  (min)  for  *ANS  file  generation 

FFP 

Boolean.  TRUE:  Force  first  pull  from  current  depth 

CDEPTFI  regardless  of  compartmental  gas  tensions  and 
prevailing  MPTTs 

FORCE  STOP 

(Boolean) 

TTIS 

Boolean  that  when  TRUE  causes  travel  time  to  each  stop 
except  the  first  air  and  first  O2  stop  in  an  ascent  to  be 
included  in  the  computed  stop  time 

RNTMODE 

residual  nitrogen  time  (RNT)  computation  mode  for 
production  of  Residual  Nitrogen  Time  Table 

•  RNTMODE  =  0:  The  residual  nitrogen  time  for  each 
RGD  is  computed  directly  as  the  time  required  at 
depth  to  reach  the  lower  bound  of  the  RGD  from  air 
saturation  at  surface. 

•  RNTMODE  =  1 :  The  residual  nitrogen  time  for  each 
RGD  is  calculated  as  the  sum  of  the  descent  time 
from  air  saturation  at  surface,  time  at  bottom,  and 
ascent  time  to  the  dropout  depth  (corresponding  to 
the  LST_DOMode  setting)  required  to  reach  the 
lower  bound  of  the  RGD  at  the  dropout  depth  during 
ascent. 

02TIME 

O2  breathing  time  per  cycle  in  O2/AIR  cycling  (min) 

AIRTIME 

Maximum  air  breathing  time  per  cycle  in  O2/AIR  cycling 
(min) 
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Parameter 

Description 

02TIME  F02 

Inspired  02  fraction  during  02  breathing  in  02/AIR  cycling 

CNDSDR_F02 

Conditional  SDR  F02 

When  constant  oxygen  fraction  breathing  gases  are 
processed,  CNDSDR_F02  is  the  fraction  of  oxygen  in  the 
inspired  gas  above  which  saturation/desaturation  time 
constant  ratios  (SDR  values)  other  than  1 .0  in  the  model 
parameter  input  file  are  used.  With  inspired  02  fractions 
less  than  the  CNDSDR_F02  value,  all  compartmental 

SDR  values  are  kept  at  1 .0,  regardless  of  their  settings  in 
the  model  parameter  input  file.  When  processing  constant 
oxygen  partial  pressure  breathing  gases,  compartmental 
SDR  values  specified  in  the  model  parameter  input  file  are 
always  used. 

02CEIL 

Maximum  depth  (fsw)  at  which  oxygen  breathing  can  start 
in  air  with  in-water  02  breathing  decompression 

GSWLAT 

Air-to-02  gas  switch  latency  (min) 

GS_DEAD 

Boolean.  TRUE:  Consider  gas  switch  latencies  as  gas 
exchange  dead  time 

AB_DEAD 

Boolean.  TRUE:  Consider  air  breaks  as  gas  exchange 
dead  time 

OMIT_TRVL 

Boolean.  TRUE:  Omit  travel  time  from  elapsed  0N02  and 
ONAIR  times  during  Air/02  cycling 
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After  execution  of  the  Readjni  routine,  TBLP7R  calls  the  RDIN7  routine  to  read  the 
maximum  permissible  tissue  tension  (MPTT)  table  from  the  model  parameter  input  file 
(*.mpf)  specified  in  the  first  line  of  *.pcf  input  file.  The  MPTT  table  comprises  a  matrix  of 
safe  ascent  criteria  arranged  in  rows  (one  for  each  stop  depth)  and  columns  (one  for 
each  gas  exchange  compartment).  This  convention  allows  use  of  tables  that  may  not  be 
built  from  solutions  to  simple  analytic  functions  of  stop  depth  and  requires  an  input  file 
with  separate  matrices  of  MPTT  values  for  each  supported  stop  depth  increment  (DINC) 
in  each  supported  depth  unit.  The  RDIN7  routine  opens  the  specified  model  parameter 
input  file,  then  locates  and  reads  the  MPTT  matrix  for  the  stop  depth  increment  and 
depth  units  specified  in  the  *.pcf  input  file.  The  dimensions  of  the  MPTT  matrix  in  the 
*.mpf  file  define  the  maximum  dive  depth  (number  of  rows  x  DINC)  and  number  of  gas 
exchange  compartments  (number  of  columns)  that  ensuing  algorithm  operation  will 
support.  The  RDIN7  routine  exits  with  a  FAIL  return  status,  and  the  start  procedure 
aborts  if  the  specified  *.mpf  file  is  not  found  or  the  MPTT  matrix  for  the  stop  depth 
increment  in  the  specified  depth  units  is  not  found  in  the  *.mpf  file. 

The  program  then  cycles  through  processing  of  each  profile  in  the  PCF  file  under 
control  of  the  ReadPCF,  Configure  Calculations  for  Table,  and  Profile  Generation  and 
Update  routines  schematized  in  Figures  2,  7,  and  9,  respectively. 
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*NOTE:  In  TBLP7R,  descent  and  ascent  rates  are  set  only  once  from  first 
values  encountered  in  the  first  profile  processed.  This  ensures  that  schedules 
in  the  table  produced  by  the  program  are  computed  with  the  same  rates. 


Figure  2.  ReadPCF  routine. 
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Figure  3.  Initialize  Profile  routine. 


Figure  4.  Initialize  Model  State  routine. 
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ENTER 

Initialize  OPTIONS 


SET: 

•  BTMTIM  =  FALSE 
.  NODSTP  =  FALSE 

•  LSTOP  =  FALSE 
.  OXYD  =  FALSE 

•  IW02D  =  FALSE 

.  STRT02  =  FALSE 

•  0N02  =  FALSE 

•  GSWLat  Accommodated  =  FALSE 

•  TT  =  0 

.  FRSTOP  =  FALSE 

•  RECYCLE  =  FALSE 
.  PSW02  =  FALSE 

•  DEAD_TIM  =  FALSE 
.  LASTOP  =  GLASTOP 


Figure  5.  Initialize  Options  routine. 


The  SetM  routine  (Figure  6)  configures  the  MPTT  matrix  to  implement  the  last  allowed 
decompression  stop  depth  at  the  user-specified  integer  multiple,  IDX,  of  the  stop  depth 
increment,  DINC.  In  a  Thalmann  Algorithm  MPTT  matrix,  the  MPTTs  for  the  i  =  1,  2,  .... 
ntiss  compartments  at  depth  D,  Mid,  are  offset  by  one  stop  depth  increment  (DINC)  to 
depth  D+DINC,  because  it  is  at  this  deeper  depth  that  the  M-,:d  are  used  to  assess  when 
ascent  to  depth  D  will  be  allowed.  Thus,  tor  DINC  =  10  fsw,  the  surfacing  MPTT  values 
(. Mj:o  =  MOj )  appear  in  the  10  fsw  row.  For  stops  at  depths  D  equal  to  integral  multiples  of 
DINC  (D  =  ADINC;  A  =  1,2,  3,  ...),  the  offset  MPTT  values  are  given  generally  by 


where  the  offset  values  are  designated  with  a  prime  and 

Mij3^=M°i+aiDx_  i_  (2) 

Eq.  (2)  is  Workman’s  original  expression15  for  the  MPTTs  at  depth  D  not  offset  by  DINC, 
where  at  is  a  slope  parameter  for  projecting  the  surfacing  MPPTs  to  depth. 
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The  MPTT  matrix  is  appropriately  configured  by  copying  the  surfacing  MVals  —  i.e. ,  the 
MPTTs  at  surface  —  to  the  row  in  the  MPTT  matrix  for  the  depth  of  the  last  allowed 
decompression  stop.  By  default,  the  MPTTs  at  surface  constitute  the  entries  for  the  1  X 
DINC  row  in  the  MPTT  table  to  support  a  last  allowed  decompression  stop  at  depth  = 
DINC.  The  routine  includes  provisions  when  called  in  MODE  =  0  mode  to  restore  the 
original  MPTT  matrix  as  read  from  the  *.mpf  file  to  support  last  allowed  stop  depth 
changes  during  profile  processing. 


Figure  6.  SetM  routine. 
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Figure  7.  Configure  Calculations  for  Table  routine. 
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Figure  8.  BTMTIM_Processing  routine. 
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Figure  9.  Profile  Generation  and  Update  routine. 
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The  Set_First_Stop_Depth  routine  (Figure  10)  contains  setup  and  completion 
procedures  for  the  FRSP7  routine  (Figure  1 1 ).  The  first  stop  depth,  DFS,  is  determined 
in  FRSP7  and  returned  by  Set_First_Stop  in  variable  DEPTH.  The  original  target  depth 
specified  in  the  *.pcf  instruction  is  returned  in  variable  FDEPTH,  and  the  switches  to 
enable  computation  of  stop  time  and  a  new  first  stop,  CSTIME  and  CFSTOP, 
respectively,  are  set  or  reset  as  appropriate. 


Figure  10.  Set  First  Stop  Depth  routine. 


The  FRSP7  routine,  schematized  in  Figure  1 1 ,  finds  the  depth  of  the  first  stop  (DFS)  for 
ascent  from  the  current  depth,  CDepth,  taking  the  effect  of  gas  exchange  during  ascent 
into  account.  This  is  accomplished  by  successive  tests  of  different  trial  first  stop  depths 
(MIND).  MIND  is  initially  set  equal  to  the  entering  instantaneous  ascent  depth  (IAD). 

The  IAD  is  then  recomputed  after  completion  of  a  trial  ascent  to  "MIND+DINC,"  where 
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DINC  is  the  stop  depth  increment.  If  the  new  IAD  differs  from  the  IAD  before  the  trial 
ascent,  the  trial  MIND  is  updated  to  this  new  value,  and  the  IAD  after  ascent  from 
CDepth  to  the  new  MIND  is  computed.  The  procedure  is  repeated  until  the  IAD  after 
ascent  and  MIND  are  equal. 

The  first  stop  depth  computed  by  this  procedure  is  never  deeper  than  the  entering 
CDepth  if  preceding  ascents  are  executed  with  appropriate  decompression  stops. 
However,  after  an  ascent  without  required  decompression  stops,  the  first  stop  depth 
found  in  this  routine,  the  IAD,  may  be  deeper  than  the  entering  CDepth.  In  such  cases 
the  first  stop  depth  is  arbitrarily  set  equal  to  the  entering  CDepth,  which  may  or  may  not 
be  an  integer  multiple  of  the  stop  depth  increment,  DINC.  Otherwise,  the  first  stop  depth 
returned  by  FRSP7  is  always  an  integer  multiple  of  DINC. 
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Figure  12.  Set  Next  Stop  Depth  routine. 


Figure  13.  (VAir  Cycling  routine. 
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Note  that  any  air-to-02  gas  switch  latency  is  completed  on  the  entering  GAS. 


Figure  14.  02-Air  Cycling  (2)  routine. 


Figure  15.  Compute  Stop  Time  routine. 
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The  Round  Time  UP  routine  schematized  in  Figure  16  conditions  the  decompression 
stop  time  computed  in  STIM7  (Figure  28)  to  be  the  next  largest  integer  value  for  output. 
If  travel  time  to  the  stop  is  being  included  in  the  output  stop  time  (TTIS  =  TRUE  and 
NOT  FRSTOP),  the  stop  time  returned  by  the  routine  is  the  actual  time  at  stop  for 
proper  model  state  update  in  a  subsequent  UPDT7  routine  (Figure  30)  call. 


► 


TIME  =  TIME +  TT 


Compute  actual  TIME  required  at  CDEPTH: 
.  TIME  =  TIME -TT 

Figure  16.  Round  Time  Up  routine. 


Node-by-node  model  state  data  is  stored  in  D,  T,  ZT,  GasCode,  P,  and  ICNTRLT  arrays 
in  the  Update/Record  model  state  routine  (Figure  17).  These  arrays  are  dimensioned  to 
support  storage  of  data  for  up  to  ten  different  diver  profiles  at  a  time  in  the  NDP.  The 
diver  pointer,  DspDvr,  is  not  used  in  TBLP7R,  because  only  a  single  schedule  or  diver  is 
processed  at  a  time.  The  violation  flag  is  also  used  only  in  the  NDP,  where  user-entered 
nodes  may  put  a  diver  in  violation  with  one  or  more  prevailing  tissue  tensions  in  excess 
of  an  MPTT.  Violations  do  not  occur  in  TBLP7R,  because  decompression  stops  are 
always  inserted  as  required  to  prevent  violations. 
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ENTER 
Update/Record  model  state 


UPDT7(CDEPTH,  DC,  TC,  GAS) 

Update  compartmental  gas  tensions 


Set  Violation  Flag  for 
current  end  node 


Clear  Violation  Flag 
for  current  end  node 


RCRD7 


FCNTRLT  (MODE=0) 

Determine  controlling  tissue,  ICNTRLT 


Store  end-segment  data  (row  CNTR  in  storage  array): 


DEPTH 

TIME  (Differential) 

TIME  (Elapsed) 

GAS  (GASCODE) 
Compartmental  gas  tensions 
ICNTRLT 


D([DspDvr,]  CNTR) 

T([DspDvr,]  CNTR) 

ZT([DspDvr,]  CNTR) 
GASCODE([DspDvr,]  CNTR)* 
P([DspDvr,]  CNTR,  i=1,NumTiss) 
ICNTRLTQDspDvr,]  CNTR) 


NOTE:  DspDvr  subscript 
used  only  in  NDP 


Advance  CNTR  for  next  entry: 
CNTR  =  CNTR  +  1 


*The  ANS  GASCODE  is  generated  from 
the  prevailing  NCODE  and,  depending  on 
the  CP02  setting,  P02  or  FN2  value. 


Return 


Figure  17.  Update/Record  Model  State  routine. 


The  Ascendr  routine  (Figure  18)  supports  the  NLIM7  and  FRSP7  routines  by  effecting  a 
no-stop  ascent  from  the  prevailing  depth  (CDepth)  to  the  passed  target  depth  (TDepth) 
or  the  surfacing  dropout  depth,  whichever  is  deeper,  while  accommodating  rig-specific 
depth-dependent  changes  in  inspired  gas  and  ascent  rate,  in  order  to  compute  the 
instantaneous  ascent  depth  (IAD)  after  the  ascent  in  UPDT7.  Thus,  for  ascents  to 
TDepth=0,  Ascendr  returns  with  IAD  evaluated  on  the  basis  of  compartmental  gas 
tensions  after  ascent  to  the  depth  less  than  or  equal  to  LASTOP  appropriate  to  the 
LST_DOMode  setting. 
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NOTE:  The  model  state  must  be  saved  in  the  calling  routine  before  Ascendr  is  called 
and  restored  after  the  Ascendr  call.  This  model  state  save/restore  is  not  handled  in 
Ascendr  because,  in  principle,  the  model  state  includes  the  computed  IAD  for  return  to 
the  calling  routine. 
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Figure  18.  Ascendr  routine. 
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As  indicated  in  Figure  18,  the  Ascendr  routine  is  supported  in  the  NDP  by  a  UBA  Gas 
Transition  Handler  routine  (Figure  19),  which  automatically  imposes  depth-dependent 
changes  in  gas  composition  appropriate  to  the  specific  “gas  mix”  or  NGAS  on  which  the 
diver  is  ascending.  Because  this  function  is  not  supported  in  TBLP7R,  depth-dependent 
changes  in  gas  composition  must  be  specified  explicitly  by  records  in  the  *.pcf  input  file. 
As  a  result,  first  decompression  stop  depths  determined  by  the  FRSP7  routine  and  no¬ 
stop  limits  determined  by  the  NLIM7  routine  in  TBLP7R  are  correct  only  if  no  depth- 
dependent  change  in  gas  composition  occurs  with  the  gas  mix  before  arrival  at  the 
dropout  depth  for  surfacing  —  i.e.,  only  if  the  ascent  transition  depth  (ATD)  is  at  or 
shallower  than  the  dropout  depth  for  surfacing. 


Figure  19.  UBA  Gas  Transition  Handler  routine  (NDP  only). 
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Depending  on  context,  three  separate  routines  handle  depth-dependent  ascent  rates  by 
parsing  ascents  into  series  of  subascents  in  depth  increments  with  ascent  rates  as 
specified  in  the  DDARATE  matrix  loaded  from  the  DDAscent.dat  file.  The  DDAscent 
Handler  routine  (Figure  20)  operates  in  conjunction  with  the  DDAscent  Handler  (3) 
routine  (Figure  21 )  to  support  ascents  with  depth-dependent  ascent  rates  in  the  main 
Profile  Generation  and  Update  routine.  An  ascent  is  parsed  into  an  appropriate  series  of 
stepwise  ascents  even  with  entries  in  DDAscent.DAT  at  higher  depth  resolution  than 
DINC. 


Figure  20.  DDAscent  Handler  routine. 
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Figure  21.  DDAscent  Handler  (3)  routine. 


The  DDAscent  Handler  (2)  routine  (Figure  22)  incorporates  all  required  procedures  to 
handle  depth-dependent  ascent  rates  in  the  Ascendr  routine. 
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ENTER  ' 
vDDAscent  Handler  (2)  / 


NOTE:  DDARATE  is  an  mx2  array  in  which 
DDARATE(m,2)  is  the  RATE  for  ascents  from 
DEPTHS  DDARATE(m,1)  to  DDARATE(m-l.l) 


Return  with  NO 
RATE  and 
entering  TDEPTH 
unchanged 


DO  m  =  1  to  100,  STEP  1 


"DDARATE(m,1)^0 


CDEPTH  <  DDARATE(m,1) 


RATE  =  DDARATE(m,2) 


DDARATE(m,2)  ^  0 


RATE  =  DDARATE(m,2) 


YES  - 

Recycle  and  - ►  MIDX  =  MIDX-1 

MIDX  >  0  „  _ _ _ _ 


TSTDEPTH  =  0 


TSTDEPTH  =  DDARATE(MIDX,1) 


TDEPTH  =  TSTDEPTH 


TSTDEPTH  >  Entering  TDEPTH 


TDEPTH  =  Entering  TDEPTH 


Figure  22.  DDAscent  Handler  (2)  routine. 


ENTER 

Output  Procedures 

TBLP7R 


Set  pointers  to  nodes 

•  end  of  deepest  isobaric  segment,  IPTRDD 

•  beginning  of  first  decompression  stop,  IPTRFS 

•  end  of  first  decompression  stop,  IPTREFS 

•  end  of  last  decompression  stop,  IPTRELS 

Depth/BT  Filter 


Skip  Profile 


NO 


NO 

> 

Total  Ascent  and 
Dive  Time  Filters 

>◄ - 

Compute  repetitive  gp  designators 

REPET7  (MODE=1,  REPETGP,  NRGD) 


YES 


'  SUBNOD  > 

Tno 

_ YES 

<RERUN^> - 

"no 

^Decompression 
,stop(s)  in  profile?, 

"  "no 


YES 


Backspace  in  PCF  file 
to  beginning  of  profile 


New  dive  depth 


NO 


YES 


Overwrite  previously  recorded  no¬ 
stop  profiles  in  the  dive  depth 

group  until  first  profile  wf  a _ 

decompression  stop.  OVRWRT  =  TRUE 


CRNDL 

NO 


YES 


Set  RERUN  = 


Record  no-stop  limit 
profile  for  the  dive 
depth  group. 


First  profile  w/  a 
decompression 
stop  in  the  dive 
depth  group. 
Compute  and 
record  no-stop  limit 
profile  for  the  dive 
depth  group. 


OVRWRT  =  FALSE 

I 


Write  Schedule  to  Main  Table  file 
[FCNTRLT(MODE=1)  for  RGD  output  suppression] 

I 


SKIPANS  : 

Tyes 


Translate  schedule  into  Augmented  NMRI 
Standard  format  &  write  to  *.ANS  file. 

Output  ANS 


NO 


Just  completed  recording  first  no-stop  limit 
profile  in  the  dive  depth  group.  Reprocess 
profile  as  originally  entered.  (Rqd 
backspace  in  PCF  file  completed  above.) 


Set  RERUN  =  FALSE  >- 


_  TRG>  NDTIME ... 


NO 


YES 

Write  Depth  Group  row  in 
SUBNOD  Table  file 

* 

REPETGP  =  0 

Backspace  in  PCF  file 

NDTIME  =  0.0 

to  beginning  of  profile 

Set  SUBNOD  =  FALSE 


Continue  SUBNOD  Table 
processing  w/hext  REPETGP 
at  current  DEPTH. 


Figure  23.  TBLP7R  Output  Procedures  routine. 


TRUE 
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The  Output  ANS  routine  (Figure  24)  supports  the  TBLP7R  Output  Procedures  routine 
and  the  Save  ANS  File  routine  in  the  NDP  by  rendering  the  profile  for  a  specified  diver 
(the  last  processed  in  TBLP7R  and  the  profile  for  passed  DspDvr  in  the  NDP)  into 
Augmented  NMRI  Standard  (ANS)  format  and  writing  the  transcribed  profile  to  an  open 
ASCII  file  for  interfacing  with  probabilistic  models  of  decompression  sickness  and  time 
of  occurrence. 


f  ENTER  N. 


Figure  24.  Output  ANS  routine. 
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Figure  25.  Configure  ANS  Data  routine. 
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Figure  26.  Set  Initial  ANS_SatGas  routine. 


The  NLIM7  routine  (Figure  27)  finds  the  remaining  no-stop  time  at  the  prevailing 
CDepth.  A  trial  no-stop  time  is  first  computed  as  the  shortest  time  required  to  saturate 
any  of  the  tissue  compartments  to  its  maximum  permissible  tissue  tension  (MPTT)  at 
storage  depth,  assuming  instantaneous  ascent  to  the  storage  depth.  This  time  is  then 
optimized  to  find  the  longest  time  that  can  be  spent  at  CDepth  before  a  no-stop  ascent 
to  storage  depth  can  be  made  at  RATE  on  GAS  without  violation.  The  No-stop  time  is 
optimized  by  using  subroutine  Ascendr  (Figure  18)  to  find  the  instantaneous  ascent 
depth  (IAD)  after  a  trial  tissue  update  at  CDepth  for  the  trial  no-stop  time  and  a 
subsequent  no-stop  ascent  to  the  storage  depth  +  DINC.  The  trial  no-stop  time  is 
increased  in  0.1 -minute  increments  until  the  IAD  after  this  ascent  increases  to  the 
storage  depth  +  DINC.  This  ensures  that  the  trial  time  is  not  too  short.  The  resultant  no¬ 
stop  trial  time  is  then  decreased  in  0.005-minute  increments  until  the  IAD  first  decreases 
to  the  storage  depth.  At  this  point  the  final  no-stop  time  is  at  most  0.005  minutes  too 
long. 
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Figure  27.  NLIM7  routine. 
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The  STIM7  routine  (Figure  28)  finds  the  time  required  at  the  current  depth,  CDepth, 
before  a  violation-free  ascent  to  Depth  can  be  executed.  Such  an  ascent  is  possible 
only  after  time  has  been  spent  at  CDepth  sufficient  to  wash  out  all  tissues  to  gas 
tensions  less  than  or  equal  to  the  maximum  permissible  tissue  tensions  in  the 
appropriate  row  IJ  of  array  M.  The  row  IJ  is  computed  from  the  target  Depth.  If  Depth  is 
not  an  integer  multiple  of  DINC,  IJ  is  computed  to  point  to  the  row  in  MPTT  array  M  for 
the  next  shallower  integer  multiple  of  the  stop  depth  increment,  DINC.  The  normal  mode 
is  for  the  next  stop  to  be  one  DINC  shallower  than  CDepth.  When  the  next  stop  is  more 
than  one  DINC  shallower  than  CDepth,  the  initially  computed  stop  time  is  shortened  in 
the  Stop  Time  Optimization  routine  (Figure  29)  to  take  advantage  of  the  compartmental 
gas  washout  that  occurs  during  ascent  from  CDepth  to  Depth  +  DINC. 
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Figure  28.  STIM7  routine. 
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Figure  29.  Stop  Time  Optimization  routine. 

The  UPDT7  routine  (Figure  30)  updates  the  tissue  inert  gas  tensions  in  array  P  for  a 
specific  depth  change  DC  over  a  specific  time  interval  TC.  Zero  (0.0)  values  for  DC  and 
TC  are  valid.  Gas  uptake  and  elimination  are  exponential  until  the  total  tissue  gas 
tension  exceeds  ambient  pressure  by  the  gas  phase  overpressure,  PBOVP,  at  which 
point  gas  elimination  becomes  time-linear.  If  PBOVP  is  set  sufficiently  high,  no 
crossover  to  time-linear  mode  occurs,  and  gas  exchange  is  always  exponential.  The 
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location  of  a  PVSAT  update  error  in  original  EL_DCA  code1,5  is  indicated  with  an 
asterisk.  The  erroneous  expression  PVSAT  =  PVSAT  +  RINRT*TNODE  has  been 
replaced  with  the  correct  expression,  PVSAT  =  PVSAT  +  RAMB  *  TNODE  . 

ASYM  models,  in  which  different  exponential  time  constants  apply  for  gas  uptake  and 
gas  elimination,  are  also  supported  in  this  routine  through  nonunitary  values  of 
SDR(TISS)  [SDR(TISS)  t  1 .0],  The  transition  between  the  two  time  constants  is  always 
made  at  a  compartmental  gas  tension  minimum  or  maximum  to  avoid  introducing  any 
discontinuity  in  the  slope  of  the  exponential  function. 
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Figure  30.  UPDT7  routine. 
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Figure  31.  UPDT7  Initialize  routine. 
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Figure  32.  Set  Time  Constants  routine. 
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The  Exponential-Linear  Crossover  Time  Computation  Procedure  (Figure  33)  computes 
TXOVER,  the  time  at  which  the  tissue  being  processed  will  transition  from  the  initial 
exponential  desaturation  mode  to  the  linear  desaturation  mode.  This  crossover  occurs 
at  the  time  when  the  total  tissue  gas  tension  exceeds  the  total  ambient  pressure  by  the 
gas  phase  overpressure  (PBOVP).  TXOVER  is  initially  assumed  equal  to  the  entering  Tl 
set  in  the  calling  routine,  UPDT7.  In  constant  P02  mode  (CP02  =  TRUE),  the  rate  of 
ambient  pressure  change,  RAMB,  and  the  rate  of  inert  gas  tension  change,  RINRT,  are 
equal,  and  TXOVER  can  be  found  analytically.  In  constant  inert  fraction  mode  (CP02  = 
FALSE),  RAMB  and  RINRT  are  not  equal,  and  TXOVER  must  be  found  by  numerical 
iteration.  In  the  latter  case,  the  constant  P02  solution  is  used  as  a  starting  value  for  the 
iterative  solution. 


Figure  33.  Exponential-Linear  Crossover  Time  Computation  routine. 
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The  Linear-Exponential  Crossover  Time  Computation  routine  (Figure  34)  computes 
TXOVER,  the  time  at  which  the  tissue  being  processed  will  transition  from  linear 
desaturation  mode  to  exponential  desaturation  mode.  This  crossover  occurs  when  the 
total  dissolved  inert  gas  tension  falls  to  a  level  that  just  exceeds  the  ambient  pressure 
by  the  gas  phase  overpressure,  PBOVP.  One  of  two  solutions  for  TXOVER  is 
applicable,  depending  on  whether  the  diver  is  in  constant  P02  mode  or  constant  inert 
fraction  mode. 


Figure  34.  Linear-Exponential  Crossover  Time  Computation  routine. 
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Figure  35.  Linear  Update  routine. 


Figure  36.  Exponential  Update  routine. 
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The  Update  IAD  routine  (Figure  37)  is  called  by  the  UPDT7  routine  to  determine  the 
new  instantaneous  ascent  depth  (IAD)  after  the  model  state  has  been  updated  to  a  new 
profile  node.  The  IAD  is  the  deepest  integer  multiple  of  the  stop  depth  increment,  DINC, 
at  which  the  gas  tension  in  any  tissue  exceeds  its  maximum  permissible  tissue  tension 
(MPTT).  As  its  name  implies,  the  IAD  is  the  shallowest  depth  to  which  an  instantaneous 
ascent  can  be  made  without  stopping  at  an  intermediate  depth  to  avoid  an  ascent 
violation.  The  IAD  is  the  surface  if  no  tissue  gas  tension  exceeds  its  shallowest  MPTT. 


MVAL  matrix,  M(i,j),  rows  assumed 
ordered  in  increasing  depth  with  j 


Figure  37.  Update  IAD  routine. 
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ENTER 

FCNTRLT  (MODE) 


ICNTRLT  =  0 


ICNTRLT  is  slowest  tissue 
with  prevailing  gas  tension 
in  excess  of  its  MVAL. 


ij  =  (CDEPTH/DINC)+1 
IDX  =  NINT(LASTOP/DINC) 
IDX=MAX(ij,IDX) 


DO  i  =  1.NTISS 


ICNTRLT  is 
slowest  tissue 
with  gas  tension 
in  excess  of  its 
MVAL  at 
beginning  of  last 
stop. 


Finding  ICNTRLT  for 
RGD  output  suppression 
in  TBLP7R 


NO 


RETURN 


DO  i  =  CNTR  -  1 ,  1,-1 


DEPTH(i)  <  LAS  I  OP 
YES 


ICNTRLT  =  Stored  ICNTRLT(i) 


END  DO 


ICNTRLT  is  tissue  with 
gas  tension  closest  to  its 
MVAL  at  end  of  last  stop. 


IDX  =  NINT(LASTOP/DINC) 


DPMIN  =oo 


DO  i  =  1.NTISS 


TDP  =  |P(i)  @  end  last  stop  -  M(i,IDX)| 


DPMIN  =TDP 
ICNTRLT  =  i 


-c 


END  DO 


NOTE:  M(i,j)  assumed  in 
order  of  increasing  half  time 
from  i  =  1  to  NTISS 


Figure  38.  FCNTRLT  routine. 


Subroutine  FCNTRLT  (Figure  38)  determines  the  tissue,  ICNTRLT,  that  controls  the  time 
at  a  decompression  stop.  The  appropriate  definition  of  ICNTRLT  depends  on  the 
application  specified  by  the  MODE  setting  in  the  FCNTRLT  call.  In  MODE  =  0  mode, 
ICNTRLT  is  determined  as  the  tissue  compartment  that  controls  the  time  at  a 
decompression  stop  on  the  basis  of  prevailing  compartmental  gas  tensions  at  the  start  of 
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the  stop.  In  MODE  =  1  mode,  ICNTRLT  is  determined  as  the  tissue  compartment  that 
controls  whether  a  repetitive  group  designator  (RGD)  will  be  included  with  a  computed 
decompression  schedule  in  an  output  table. 

Subroutine  FCNTRLT  is  called  in  MODE  =  1  mode  only  by  the  TBLP7R  Output 
Procedures  routine  (Figure  23)  after  each  computed  schedule  is  complete  with  the  diver 
at  surface.  In  this  mode  FCNTRLT  determines  ICNTRLT  as  the  controlling  tissue 
compartment  for  surfacing,  depending  on  the  SRF_CNTRLT_MODE  switch  setting  from 
the  initialization  file: 

•  SRF_CNTRLT_MODE  =  0:  The  controlling  compartment  for  surfacing  is  the 
highest  half  time  (slowest)  compartment  that  had  a  gas  tension  in  excess  of  its 
surfacing  M-Val  when  the  diver  arrived  at  the  last  decompression  stop  depth 
during  ascent. 

•  SRF_CNTRLT_MODE  =  1 :  The  controlling  compartment  for  surfacing  is  the 
compartment  that  had  gas  tension  closest  to  its  surfacing  M-Val  when  the  diver 
left  the  last  decompression  stop  depth  during  ascent. 

An  RGD  is  then  included  with  the  schedule  in  the  output  table  only  if  ICNTRLT  is  less 
than  or  equal  to  the  reference  tissue,  IREF  (i.e.,  only  if  ICNTRLT  points  to  a  tissue 
compartment  with  a  gas  exchange  half-time  faster  than  or  equal  to  that  of  the  reference 
tissue). 

Subroutine  FCNTRLT  is  called  in  MODE  =  0  mode  in  the  Update/Record  Model  State 
routine  (Figure  17)  in  TBLP7R,  or  at  the  end  of  the  STIM7  routine  (Figure  28)  in  the 
NDP,  as  each  decompression  stop  is  processed.  In  the  NDP,  the  returned  value  of  the 
controlling  tissue,  ICNTRLT,  is  then  used  in  the  REPET7  routine  (Figure  39)  to 
determine  the  prevailing  RGD  for  display.  A  value  of  ICNTRLT  =  0  is  returned  for  no¬ 
stop  ascents  reflecting  that  no  tissue  controls  because  all  tissue  gas  tensions  during 
traverse  through  the  last  stop  depth  are  less  than  corresponding  surfacing  MVALs). 


The  RGD  for  each  schedule  is  determined  in  subroutine  REPET7  (Figure  39).  Original 
EL-DCA  subroutines  for  computation  of  surfacing  repetitive  groups,  surface  interval 
credit  tables,  and  residual  nitrogen  times  for  repetitive  diving  are  described  in 
Thalmann.4  Detailed  code  for  the  modules  from  which  the  REPET7  routine  is  adapted  is 
given  in  an  accompanying  untitled  technical  memorandum.5  The  diver  is  presumed  to 
be  at  surface  when  this  routine  is  called.  REPET7  is  only  called  in  MODE  =  0  mode  in 
the  NDP. 
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Figure  40.  Compute  RGT  routine. 
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ENTER 

Ascnd4RGT 


LASTM  =  FALSE 


""  LASTM  & 
PREF  <  PFINL 


PREF  >  PFINL 

T  YES 


TRG  =  TRG  +  0.1 


TRG  =  MAX(  TRG  -  0.005,0  ) 


NO  ,  '  \  YES 

-  TRG  >  0  Y - 


LASTM  =  TRUE 


Figure  41.  Ascnd4RGT  routine. 


3.  Composite  Table  Assembly  (BLDCMPTBL) 

The  Air  Decompression  Table  in  the  U.S.  Navy  Diving  Manual,  Revision  6,  is  a  single 
table  comprising  schedules  for  three  modes  of  decompression  for  each  tabulated  dive 
depth  and  time  combination:  (1)  in-water  decompression  on  air,  (2)  in-water 
decompression  on  air  and  oxygen,  and  (3)  surface  decompression  with  oxygen 
(SurD02).  For  each  tabulated  dive  depth  and  time  combination,  the  in-water  oxygen 
decompression  schedule  appears  in  a  line  beneath  the  in-water  air  schedule,  with  stop 
times  that  indicate  only  the  oxygen  breathing  time  after  arrival  at  the  02CEIL  depth  of 
30  fsw  and  a  total  ascent  time  that  includes  all  air  breathing  time  during  any  air- 
breathing  breaks  scheduled  after  start  of  in-water  oxygen  breathing.  Each  pair  of  in¬ 
water  air/in-water  oxygen  schedules  is  accompanied  by  a  repetitive  group  designator 
(RGD)  for  surfacing  after  either  type  of  decompression  (taken  as  the  highest  RGD 
computed  for  the  pair)  and  an  entry  that  gives  the  number  of  oxygen  breathing  periods 
in  the  SurD02  schedule  for  the  dive.  The  latter  is  computed  by  rule  from  the  oxygen 
breathing  time  in  the  in-water  oxygen  decompression  schedule.11 

The  BLDCMPTBL  program  (Figure  42)  was  developed  to  compile  the  composite  table 
and  the  accompanying  ANS-coded  files  for  the  tabulated  schedules  from  the  output  files 
of  a  TBLP7R  run  on  each  of  three  different  PCF  input  files.  The  input  files  must  contain 
profile  templates  for  air  dives  to  the  same  depths  and  bottom  times  but  specify  different 
options  for  completing  the  schedules: 

PCF  file  (1)  -  no  options:  yields  schedules  for  in-water  decompressions  on  air. 

PCF  file  (2)  -  GIW  global  option  and  AB_DEAD  parameter  set  TRUE:  yields 

schedules  for  in-water  oxygen  decompressions  with  air-breathing 
breaks  during  oxygen  breathing  periods. 

PCF  file  (3)  -  GIW  and  NAB  global  options:  yields  schedules  for  in-water  oxygen 

decompressions  with  no  air-breathing  breaks  during  oxygen  breathing 
periods. 

With  the  AB_DEAD  parameter  set  TRUE,  in-water  oxygen  breathing  times  in  schedules 
produced  from  PCF  file  (2)  are  the  same  as  the  in-water  oxygen  breathing  times  in 
schedules  produced  from  PCF  file  (3),  but  only  the  schedules  produced  from  PCF  file 
(2)  list  stop  times  that  include  air  breathing  times  in  air-  breaks  during  in-water  oxygen 
decompression.  Schedules  produced  from  PCF  file  (2)  consequently  list  the  correct  total 
ascent  time  for  each  profile,  but  schedules  produced  from  PCF  file  (3)  list  the  in-water 
oxygen  breathing  times  required  to  compute  corresponding  SurD02  schedules. 


Parameters  specific  to  the  calculation  of  SurD02  schedules  (Table  2)  are  input  by  the 
Readjni  routine  from  a  [CMPTBL  SURD02]  section  of  the  DCA_Params.ini  file  during 
BLDCMPTBL  startup  (Figure  42). 
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Figure  42.  BLDCMPTBL  Main  routine. 
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Table  2.  SurD02  Parameters 

(Loaded  by  READJNI  routine  from  [CMPTBL  SURD02] 
section  of  program  initialization  file  at  program  startup] 


Parameter 

Description 

AIRTIME 

Maximum  air  breathing  time  (min)  in  an  oxygen-air  cycle  during 
oxygen  breathing  in  SurD 

CARATE 

Chamber  ascent  rate 

CDRATE 

Chamber  descent  rate 

DrpOut_ARATE 

Ascent  rate  for  dropout  to  surface  during  SurD  initiation 

DrpOut_DEPTH 

Depth  from  which  to  initiate  SurD  and  dropout  to  surface 

SurDTimFctr 

Multiplier  for  computed  O2  breathing  times  in  SurD02  schedules 

02TIME 

Maximum  oxygen  breathing  time  (min)  in  an  oxygen-air  cycle 
during  oxygen  breathing  in  SurD 

02TIME_F02 

Fraction  of  oxygen  in  inspired  gas  during  oxygen  breathing  in 

SurD 
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ENTER 

File  Read  and  Synch 


H  Read  TBL1  record:  INSTR1 


INSTR1  is  a  schedule  line 


PASS  =  FALSE  "> 


INSTR2  is  a  schedule  line 


H  Read  TBL.3  record:  INSTR3  U 


INSTR3  is  a  schedule  line 


READ  schedule  dive  depth  D1  and  bottom  time  BT1  from  INSTR1 
READ  schedule  dive  depth  D2and  bottom  time  BT2from  INSTR2  « 
READ  schedule  dive  depth  D3and  bottom  time  BT3from  INSTR3 


►<:"  PASS2  =  FALSE 


PASS3=  FALSE 


BT2  =  BT1_cT> 


PASS2  =  TRUE  > 


If  NOT  PASS,  COPY  ANSI  profile  into 
ANSI  buffer  (Overwrite  any  preexisting) 


COPY  ANS2  profile  into  ANS2  buffer 
(Overwrite  any  preexisting) 


PASS  =  TRUE 


PASS3  =  TRUE 


D2  >  D3  ;> 


Figure  43.  File  Read  and  Synch  routine. 


Edit  ISTR3 


Edit  ISTR1 


Figure  44.  ASSMBL  CMPTBL  routine. 
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Figure  45.  Gen  SurD02-ANS  routine. 
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Figure  46.  Chamber  Decompression  routine. 
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4.  Surface  Interval  Credit/Residual  Nitrogen  Time  Table  Generator  (SICT7) 

The  SICT7  program  was  first  described  in  conjunction  with  development  of  new 
repetitive  and  multilevel  dive  procedures  and  tables  for  MK  16  MOD  0  0.7  ata  constant 
P02-in-N2  diving.4  FORTRAN  code  for  the  original  version  of  SICT7  was  published  in 
NEDU  TR  09-85. 5  In  the  new  procedures,  residual  nitrogen  times  (RNTs)  were 
determined  from  the  Sub-NoD  and  Main  Decompression  Tables  via  end-interval  bottom 
number  groups,  which  obviated  need  for  accompanying  RNT  tables.  The  original  SICT7 
was  consequently  designed  to  generate  only  surface  and  shallow  interval  credit  tables 
with  end-interval  bottom  number  group  repetitive  group  designators.  The  bottom 
number  group  method  for  determining  RNTs  was  not  adopted,  however.  Reversion  to 
the  letter  group  method  for  determining  RNTs  for  MK  16  MOD  0  0.7  ata  constant  P02- 
in-N2  diving  required  conversion  of  final  number  repetitive  group  designators  to  their 
approximate  letter  values  in  the  surface  interval  credit  table  and  generation  of  an 
accompanying  RNT  table.  Lacking  the  software  to  generate  the  latter  table  directly,  the 
table  was  assembled  from  RNTs  in  the  Sub-NoD  Table,  each  of  which  includes  descent 
time  to  bottom  depth  and  the  influence  of  gas  washout  during  ascent  to  the  last  allowed 
decompression  stop  depth  after  completion  of  the  bottom  time.  The  capability  to 
compute  surface  and  shallow  interval  offgassing  credits  for  letter  group  repetitive 
designators  and  to  generate  an  RNT  table  directly  for  the  surface  interval  credit  table 
was  incorporated  into  SICT7  only  later  for  the  development  of  tables  for  MK  16  MOD  1 
1 .3  ata  constant  P02-in-N2  diving6  and  MK  16  MOD  1  1 .3  ata  constant  P02-in-He 
diving.7,8  The  enhanced  SICT7  generates  an  RNT  table  with  entries  for  each  repetitive 
group  equal  to  the  corresponding  entries  in  the  Sub-NoD  Table  with  RNTMODE  =  1 . 
With  RNTMODE  =  0,  the  enhanced  SICT7  generates  an  RNT  table  with  entries  for  each 
repetitive  group  that  do  not  include  descent  times  or  the  influences  of  gas  washout 
during  subsequent  ascent  (see  Compute  RGT  routine,  Figure  40). 

Similar  to  TBLP7R,  SICT7  operates  by  processing  profile  templates  read  from  a  PCF 
file.  However,  for  RNT  table  generation,  this  file  should  contain  only  one  profile  template 
for  each  repetitive  dive  depth  desired  in  the  final  RNT  table,  and  the  profile  templates 
should  be  in  order  of  increasing  dive  depth.  As  with  Sub-NoD  table  generation,  dive 
bottom  times  in  the  PCF  profiles  are  ignored. 

SICT7  operation  starts  with  execution  of  the  procedures  illustrated  in  Figure  47. 
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Figure  47.  SICT7  Main  routine. 
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The  Generate  SICT  routine  (Figure  48)  computes  a  Surface  Interval  Credit  Table  or  a 
Shallow  Interval  Credit  Table  depending  on  the  entering  value  of  IDEPTH  set  in  the 
calling  SICT7  Main  routine.  Tables  are  generated  with  starting  letter  repetitive  group 
designators  if  IDEPTH  =  1  or  2,  or  with  starting  bottom  number  repetitive  group 
designators  if  IDEPTH  =  3,  4,  or  5.  Tables  with  starting  letter  repetitive  group 
designators  may  have  either  letter  or  bottom  number  final  repetitive  group  designators 
depending  on  the  value  of  the  “Use  Letter  Group”  (ULG)  switch  set  in  accord  with  user 
response  to  a  query  at  program  start.  In  any  case,  results  for  a  given  table  are  stored  in 
a  32  x  32  RTIM(J,K)  array.  The  starting  tissue  tension  PSTART  for  each  starting  group 
is  the  highest  reference  tissue  tension  in  the  group.  The  final  tissue  tension  PFINL  for 
each  final  group  is  the  highest  reference  tissue  tension  in  the  next  lower  group.  The 
letter  repetitive  group  gas  tension  increment  RGI  (see  SICT  Setup  routine,  Figure  49) 
and  all  PSTART  and  PFINL  values  are  defined  with  respect  to  the  saturation  inert  gas 
tension  in  the  reference  compartment  at  surface,  PINF,  offset  by  addition  of  a  1  fsw 
allowed  depth  error  documented  on  page  5  of  Thalmann.4 

Starting  inert  gas  tensions  in  the  reference  compartment  for  each  starting  repetitive 
group  designator  in  the  table  (j  =  1,  NGRU)  are  stored  in  the  (j,  1 )  partition  of  the  RTIM 
array.  Final  inert  gas  tensions  in  the  reference  compartment  for  the  kth  final  repetitive 
group  designator  in  the  table  (k  =  1 , 31 )  are  stored  in  the  (32,k+1 )  partition  of  the  RTIM 
array.  The  times  at  the  prevailing  CDEPTH  required  to  reach  the  kth  final  group  gas 
tension  PFINL  from  the  jth  starting  group  tension  PSTART  are  stored  in  the  intervening 
RTIM  elements,  RTIM(j,k+1).  The  stored  final  inert  gas  tensions  and  times  are 
consequently  shifted  in  the  RTIM  array  by  one  group  to  the  right,  so  that  the  (32, k) 
element  is  the  final  inert  gas  tension  and  the  (j,k)  elements  are  the  corresponding  end- 
interval  times  for  the  next  lower  k-1  group. 

On  output  to  the  product  interval  credit  table,  interval  times  stored  in  the  RTIM  array  for 
each  starting  repetitive  group  designator)  are  written  in  reverse  order.  For  tables  with 
starting  and  final  letter  repetitive  groups,  the  number  of  final  repetitive  groups,  NFGRU, 
equals  16,  the  number  of  starting  repetitive  groups,  and  the  times  are  written  to  two 
consecutive  lines.  In  the  first  line,  the  starting  letter  repetitive  group  designator  is 
followed  by  interval  start  times  beginning  with  the  value  for  the  highest  group  Z  in 
RTIM(j,NGRU+2)  and  ending  with  the  value  for  the  lowest  group  A  in  RTIM(j,3)  after  1 
min  is  added  to  each  stored  value  and  the  result  is  rounded  to  the  nearest  integer. 
Corresponding  interval  end  times  beginning  with  the  value  for  the  highest  group  Z  in 
RTIM(j,  NGRU+1)  and  ending  with  the  value  for  the  lowest  group  A  in  RTIM(j,2)  are 
written  to  the  second  line  after  each  stored  value  is  rounded  to  the  nearest  integer. 

For  tables  with  final  bottom  number  repetitive  groups,  the  number  of  final  repetitive 
groups,  NFGRU,  is  21  for  IDEPTH  =  1 , 2,  or  5  and  31  for  IDEPTH  =  3  or  4,  and  only  a 
single  line  giving  end-interval  times  is  written  for  each  starting  repetitive  group 
designator).  The  starting  repetitive  group  designator  [letter  for  a  surfacing  letter  group 
or  stored  starting  gas  tension  in  RTIM(j,1 )  for  a  bottom  number  group]  is  followed  by 
interval  end  times  beginning  with  the  value  for  the  highest  group  in  RTIM(j,  NFGRU+1 ) 
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and  ending  with  the  value  for  the  lowest  group  in  RTIM(j,2)  after  each  stored  value  is 
rounded  to  the  nearest  integer. 

After  the  interval  times  for  all  starting  groups  have  been  written,  the  NFGRU  final  group 
designators  are  written  in  reverse  order  to  the  last  line  in  the  table.  Final  bottom  number 
group  designators  are  the  stored  final  gas  tensions  in  the  (32, k)  [k  =  3,  NFGRU+2] 
partition  of  the  RTIM  array.  The  number  group  designator  written  in  the  left-most 
column,  RTIM(32, NFGRU+2),  is  undefined  but  taken  as  the  value  of 
INT(RTIM(32,NFGRU+1  ))+1 .  Because  the  stored  final  gas  tension  in  each  (32, k) 
element  of  the  RTIM  array  is  the  final  tension  for  the  (k-1 )  final  group,  use  of  the  k  =  3, 
NFGRU+2  elements  shifts  the  final  number  group  designators  one  column  to  the  right  in 
the  final  table,  which  aligns  each  designator  with  the  column  that  contains  the  maximum 
interval  end  times  to  which  the  designator  applies. 
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Figure  48.  Generate  SICT  routine. 


62 


Figure  49.  SICT  Setup  routine. 
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The  Surface  Interval  Credit  Table  produced  by  the  Generate  SICT  routine  with  IDEPTH 
=  1  is  accompanied  by  a  Residual  Nitrogen  Time  (RNT)  Table  produced  by  the 
Generate  RNT  Table  routine  schematized  in  Figure  50. 


Figure  50.  Generate  RNT  Table  routine. 


Each  row  jj  of  residual  gas  times  in  the  RNT  Table  is  computed  in  the  Compute  Row 
routine  (Figure  51 )  with  the  jjth  profile  template  in  the  PCF  input  file  specified  by  the  user 
during  program  startup  (Figure  47).  The  repetitive  dive  depth  for  the  row  is  stored  in  the 
jjth  element  of  an  RNTD  vector,  and  the  corresponding  residual  gas  times  for  the 
NFGRU  final  repetitive  groups  are  stored  in  the  (jj,kk)  [kk=1 ,31]  partition  of  an  RNT 
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array.  For  all  jj  and  j,  the  (jj,kk)  elements  of  the  RNT  array  correspond  to  the  (j,kk+1 ) 
elements  of  the  RTIM  array  in  the  Generate  SICT  routine  (Figure  48). 

After  the  residual  gas  times  for  all  repetitive  dive  depths  in  the  PCF  input  file  have  been 
computed,  the  product  RNT  table  is  written  under  its  corresponding  interval  credit  table 
with  times  for  each  depth  aligned  in  a  row  under  the  final  repetitive  group  columns  to 
which  they  apply. 
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Figure  51.  Compute  Row  routine. 
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Note  that  the  RNT  for  each  end-interval  number  group  is  computed  as  the  time  at 
CDEPTH  to  reach  the  reference  tissue  tension  of  the  next  higher  group.  For  each  end- 
interval  letter  group,  the  RNT  is  computed  as  the  time  to  reach  the  highest  reference 
tissue  tension  in  the  group.  This  causes  each  RNT  table  entry  to  give  the  most 
conservative  or  maximum  time  penalty  for  its  corresponding  end-interval  final  group. 


Figure  52.  ReadPCF  SICT7  routine. 
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6.  Appendix  A.  List  of  Variables 
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AIRTIME 

AMBA02 

ANSTDiff 

ASTIM 

BLANK 

BTMAX 

BTMTIM 

CDEPTH 

CF 

CFSTOP 

CNDSDR_F02 

CNTR 

CONC 

CONLBL 

CP02 

CSTIME 

D 

DAYTIM 

DC 

DDARATE 

DDASCENT 

DDGDEPTH 

DEPTH 

DFS 

DINC 

DONE 

DSTOPS 

DUMP 

EOF 

FDEPTH 

FFP 

FN2 

FRSTOP 


FULBUF 

GAS 


AIR  BREATHING  TIME  (MIN)  IN  AN  02-AIR  CYCLE 
AMBIENT-ARTERIAL  OXYGEN  GRADIENT  (FSW) 

BOOLEAN.  TRUE:  DIFFERENTIAL  TIMES  IN  ANS  OUTPUT 
FALSE:  ELAPSED  TIMES  IN  ANS  OUTPUT 
TOTAL  ASCENT  TIME  (MIN) 

DUMMY  CHARACTER  PLACEHOLDER 
MAXIMUM  PERMISSIBLE  BOTTOM  TIME  (MIN) 

BOOLEAN.  TRUE:  INCLUDE  DESCENT  TIME  IN  'TIME' 
CURRENT  DEPTH  (FSW  OR  MSW) 

METERS  SEAWATER  (MSW)  TO  FEET  SEAWATER  (FSW) 
CONVERSION  FACTOR  (=  1.0/MSWpFSW) 

BOOLEAN.  TRUE:  CALCULATE  FIRST  STOP 
THRESHOLD  FI02  FOR  USE  OF  SDR(I).NE.l 
POINTER  TO  MODEL  STATE  "TT"  ARRAY  FOR  UPDATES 
IN  UPDATE/RECORD  MODEL  STATE  SUBROUTINE 
INSPIRED  OXYGEN  TENSION  FOR  PRINTOUT 
INSPIRED  OXYGEN  TENSION  LABELS  FOR  PRINTOUT 
BOOLEAN.  TRUE:  CONSTANT  INSPIRED  02  PARTIAL 
PRESSURE;  FALSE:  CONSTANT  INSPIRED  02  FRACTION 
BOOLEAN.  TRUE:  COMPUTE  DECOMPRESSION  STOP  TIME 
PROFILE  DEPTH  ARRAY 
DATE  TIME  ARRAY 

DEPTH  CHANGE  BETWEEN  PROFILE  SUB-SEGMENT  END- 
NODES  (FSW  OR  MSW) 

DEPTH-DEPENDENT  ASCENT  RATE  DATA  ARRAY 
BOOLEAN.  TRUE:  USE  DEPTH-DEPENDENT  ASCENT  RATE 
DEPTH  OF  CURRENT  DIVE  DEPTH  GROUP  (FSW  OR  MSW) 
NEXT  DEPTH  (FSW  OR  MSW) 

DEPTH  OF  FIRST  STOP  (FSW  OR  MSW) 

STOP  DEPTH  INCREMENT  (FSW  OR  MSW) 

BOOLEAN.  TRUE:  FINISHED  ENTERING  CURRENT 
PROFILE 

DEPTHS  OF  DECOMPRESSION  STOPS  (FSW  OR  MSW) 

DUMMY  PLACEHOLDER 

END  OF  FILE  ENCOUNTERED  DURING  SEQUENTIAL  READ 
FINAL  DEPTH  SPECIFIED  IN  PCF  INSTRUCTION  (FSW 
OR  MSW) 

BOOLEAN.  TRUE:  FORCE  FIRST  PULL  TO  NEXT  STOP 
CURRENT  INERT  GAS  FRACTION  IN  USE 
BOOLEAN.  TRUE:  COMPUTING  PROPERTIES  FOR  FIRST 
STOP  IN  AN  ASCENT  OR  FOR  FIRST  STOP  AFTER  A  GAS 
SWITCH 

BOOLEAN.  TRUE:  'I PRO'  BUFFER  FULL 
PROFILE  INERT  GAS  TENSION  ARRAY 
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GASCODE 

GASCODEu 

GASTSN 

HLFTM 

IAD 

ICHNG 

ICODE 

ICONC 

I  GAS 

IGRU 

ILBL 

I  PRO 

IPRT 

IARATE 

IDRATE 

IPTRDD 

IPTRFS 

IPTREFS 

IPTRELS 

IREF 

IRGD 

ISTOP 

ITIS 

K 

KDSAT 

KSAT 

LASTOP 

LB 

LIMLIN 

LLINE 

LP 

LSTOP 

LU 

M 

MAXNGR 

MAXPRO 


AUGMENTED  NMRI  FORMAT  GASCODE  ARRAY 
AUGMENTED  NMRI  FORMAT  GASCODE  ARRAY  (SHIFTED 
FOR  OUTPUT  TO  ANS  FILE) 

GAS  TENSION  VALUE  ARRAY 
COMPARTMENT  HALFTIMES  ARRAY  (MIN) 

INSTANTANEOUS  ASCENT  DEPTH  (FSW  OR  MSW) 

BOTTOM  DEPTH  OF  PROFILE  GROUP  BEING  PRINTED 
INTEGER  PORTION  OF  CURRENT  GASCODE 
POINTS  TO  P02  OR  F02  LABELS  IN  ' CONLBL ' 

INERT  GAS(ES)  NAME (S) 

REPETITIVE  GROUP  DESIGNATOR 

POINTS  TO  UNITS  OR  SPECIFICATION  IN  'CONLBL' 
DECOMPRESSION  PROFILE  ARRAY 
PRINT  MODEL  PARAMETER  INPUT  FILE? 

FIRST  SPECIFIED  ASCENT  RATE  IN  FIRST  PROFILE 
(FSW  OR  MSW/MIN) 

FIRST  SPECIFIED  DESCENT  RATE  IN  FIRST  PROFILE 
(FSW  OR  MSW/MIN) 

POINTER  TO  ENTRY  IN  D  AND  T  ARRAYS  FOR  LAST 
NODE  AT  MAX  DIVE  DEPTH 

POINTER  TO  ENTRY  IN  D  AND  T  ARRAYS  FOR  FIRST 
NODE  AT  FIRST  DECO  STOP 

POINTER  TO  ENTRY  IN  D  AND  T  ARRAYS  FOR  LAST 
NODE  AT  FIRST  DECO  STOP 

POINTER  TO  ENTRY  IN  D  AND  T  ARRAYS  FOR  LAST 
NODE  AT  LAST  STOP 

REFERENCE  COMPARTMENT  NUMBER  FOR  REPETITIVE  GPS 
NUMERICAL  REPETITIVE  GROUP  DESIGNATOR  ARRAY 
STOP  DEPTH  AT  CURRENT  'IPRO'  POSITION  (FSW  OR 
MSW) 

LOOP  CONTROL;  COMPARTMENT  SUBSCRIPT 
CURRENT  POSITION  IN  DIVE  PROFILE  ARRAYS 
TIME  CONSTANT  FOR  DESATURATING  TISSUE  (1/MIN) 
TIME  CONSTANT  FOR  SATURATING  TISSUES  (1/MIN) 
LAST  (OR  SHALLOWEST)  ALLOWED  STOP 
PROFILE  COORDINATE  INPUT  FILE  DEVICE  NUMBER 
BOOLEAN.  TRUE:  PRINT  LIMIT  LINE  BEFORE  THIS 
SCHEDULE 

ARRAY  HOLDING  VALUE  'LIMLIN'  FOR  EACH  SCHEDULE 
DEVICE  NUMBER  FOR  LINE  PRINTER 

BOOLEAN.  TRUE:  LAST  ENTERED  DEPTH  IS  DEPTH  OF 

LAST  STOP  BEFORE  SURFACING 

DEVICE  NUMBER  FOR  TERMINAL 

MAXIMUM  TISSUE  GAS  TENSION  ARRAY  (FSW) 

MAXIMUM  NUMERICAL  REPET  GROUP  VALUE 
MAXIMUM  NUMBER  OF  PROFILES  PER  PAGE. 
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MAXSTP 

METRIC 

MINSTP 

MORE 

MPIF 

MSWpFSW 

NCODE 

NCOLUMN 

NDFORM 

NDTIME 

NDTIMEP 

NGAS 

NFGRU 

NGRU 

NODLIM 

NODSTP 

NOLGRU 

NONGRU 

NORCRD 

NOPRNT 

NPRO 

NRGD 

NRIP 

NS1-NS5 

NSPC 

NSPCD 

NTISS 

OPTN 

020NCLS 


02STRTD 

02TIME 

02TIME_F02 

OXYD 


MAXIMUM  NUMBER  OF  STOPS  ALLOWED  BY  PRINTOUT 
FORMAT 

BOOLEAN.  TRUE:  DEPTH  AND  RATE  INPUTS  IN  METERS 
SEAWATER  (MSW) 

MINIMUM  NUMBER  STOPS  PER  PAGE 

BOOLEAN.  TRUE:  PROCESS  ANOTHER  DIVE  PROFILE 

MODEL  PARAMETER  INPUT  FILENAME 

FEET  SEAWATER  (FSW)  TO  METERS  SEAWATER  (MSW) 

CONVERSION  FACTOR 

INTEGER  PORTION  OF  ANS  GASCODE 

MAXIMUM  NUMBER  OF  PRINTER  COLUMNS 

BOOLEAN.  TRUE:  PRINTOUT  IN  SUB  NO-D  FORMAT 

NO-DECOMPRESSION  TIME,  LESS  DESCENT  TIME 

NO-DECOMPRESSION  TIME,  INCLUDING  DESCENT  TIME 

GAS  TENSION  ARRAY,  ' GASTSN ' ,  SUBSCRIPT 

NUMBER  OF  FINAL  REPETITIVE  GROUPS  IN  SURFACE 

INTERVAL  CREDIT  TABLE 

NUMBER  OF  STARTING  REPETITIVE  GROUPS  IN  SURFACE 
INTERVAL  CREDIT  TABLE 

BOOLEAN.  TRUE:  COMPUTE  NO-DECOMPRESSION  LIMIT 
BOOLEAN.  TRUE:  DON'T  COMPUTE  DECOMPRESSION 
STOPS  (FORCED  NO-STOP  ASCENT) 

BOOLEAN.  TRUE:  DON'T  PRINT  LETTER  GROUP 
BOOLEAN.  TRUE:  DON'T  PRINT  NUMERICAL  GROUP 
BOOLEAN.  TRUE:  DON'T  RECORD  DIVE  PROFILE 
DUMMY  CHARACTER  READ  TO  MAINTAIN  COMPATIBILTY 
W/DMDB7 

NUMBER  OF  PROFILES  RECORDED  IN  'IPRO'. 

NUMERICAL  REPETITIVE  GROUP  DESIGNATOR 
NUMBER  OF  RECORDS  IN  PROFILE 

NUMBER  VARIABLES  FOR  VARIABLE  FORMAT  STATEMENTS 
NUMBER  OF  SPACES  TO  FIRST  STOP  FOR  A  PROFILE 
NUMBER  OF  SPACES  TO  ACCOMODATE  LONGEST  PROFILE 
NUMBER  OF  GAS  EXCHANGE  COMPARTMENTS  IN  ALGORTHM 
( 2  0  MAX . ) 

OPTION  ARRAY  FOR  CURRENT  PCF  RECORD 
BOOLEAN.  TRUE:  START  02  BREATHING  ON  CHAMBER 
LEAVE  SURFACE  DURING  SURD02 

FALSE:  START  02  BREATHING  ON  ARRIVAL  AT 
DEPTH  DURING  SURD02 

BOOLEAN.  TRUE:  OXYD=TRUE  AND  02-Air  CYCLING 
STARTED 

OXYGEN  BREATHING  TIME  (MIN)  IN  AN  02-AIR  CYCLE 
F02  DURING  02TIME  IN  02-AIR  CYCLING 
BOOLEAN.  TRUE:  USE  OXYGEN  DECOMPRESSION  W /  02- 
AIR  CYCLING 
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C  0N02 

C  PHEAD 

C  PLLINE 

C  P 

C  PAC02 

C  PAMB 

C  PAN  2 

C  PA02 

C  PBOVP 

C  PH20 

C  PINF 

C  PLIN 

C  P02 

C  PREF 


C  PROFL1 

C  PRO ID 

C  PSEP 

C  PSTART 

C  PSTOPS 

C  PTITLE 

C  PVC02 

C  PVN2 

C  PV02 

C  PVSAT 

C  RAMB 

C  RATE 


C  RECYCLE 

C  RGD 

C  RGDO 

C  RGI 

C  RGRU 

C  RGT 

C  RINRT 


BOOLEAN.  BREATHING  02  (.TRUE.)  OR  BREATHING  AIR 
(.FALSE.)  IN  02-AIR  CYCLING 

CHARACTER  VARIABLE  HOLDING  PAGE  HEADER  FORMAT 
CHARACTER  VARIABLE  HOLDING  LIMIT  LINE  FORMAT 
COMPARTMENT  INERT  GAS  TENSION  ARRAY  (FSW) 
ARTERIAL  C02  TENSION  (FSW) 

AMBIENT  HYDROSTATIC  PRESSURE  (FSW) 

ARTERIAL  INERT  GAS  TENSION  (FSW) 

ARTERIAL  02  TENSION  (FSW) 

GAS  PHASE  OVERPRESSURE  (FSW) 

PARTIAL  PRESSURE  OF  WATER  VAPOR  (FSW) 

TISSUE  INERT  GAS  TENSION  WHEN  SATURATED  ON  AIR 
AT  SURFACE 

COMPARTMENTAL  INERT  GAS  TENSION  AT  TXOVER  AND 
DURING  LINEAR  UPDATE  (FSW) 

PREVAILING  INSPIRED  OXYGEN  PARTIAL  PRESSURE 
(ATA) 

INERT  GAS  TENSION  IN  REFERENCE  COMPARTMENT  AT 
END  SHALLOWEST  STOP  FOR  LETTER  REPETITIVE  GROUPS 
OR  AT  END  BOTTOM  TIME  FOR  NUMBER  REPETITIVE 
GROUPS 

BOOLEAN.  TRUE:  FIRST  PROFILE  TO  BE  READ 
PROFILE  IDENTIFICATION  STRING 

CHARACTER  VARIABLE  FOR  PROFILE  SEPARATOR  FORMAT 
STARTING  INERT  GAS  TENSION  IN  REFERENCE 
COMPARTMENT  FOR  REPETITIVE  GROUP  CALCULATION 
CHARACTER  VARIABLE  HOLDING  PROFILE  FORMAT 
CHARACTER  VARIABLE  HOLDING  HEADER  FORMAT 
VENOUS  C02  TENSION  (FSW) 

VENOUS  INERT  GAS  TENSION  (FSW) 

VENOUS  02  TENSION,  PREVAILING  (FSW) 

VENOUS  INERT  GAS  TENSION  WHEN  SATURATED  AT 
PREVAILING  PAMB  (FSW) 

RATE  OF  AMBIENT  PRESSURE  CHANGE  IN  CURRENT 
PROFILE  SUB-SEGMENT  (CONVERTED:  FSW/MIN) 

RATE  OF  AMBIENT  PRESSURE  CHANGE  IN  CURRENT 
PROFILE  SUB-SEGMENT  (UNITS  IN  PCF  FILE:  FSW/MIN 
OR  MSW/MIN) 

BOOLEAN.  PROCESS  FLOW  CONTROL 

NUMERICAL  REPETITIVE  GROUP  DESIGNATOR  =  NRGD 
LETTER  REPETITIVE  GROUP  DESIGNATOR 
INERT  GAS  TENSION  INCREMENT  BETWEEN  SUCCESSIVE 
LETTER  REPETITIVE  GROUPS 

LETTER  REPETITIVE  GROUP  DESIGNATOR  ARRAY 
TIME  TO  REACH  REPETITIVE  GROUP 
RATE  OF  INSPIRED  INERT  GAS  TENSION  CHANGE  IN 
CURRENT  PROFILE  SUB-SEGMENT  (FSW/MIN) 


A-4 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 
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c 
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R02 

RTIM 

RTIME 

SDR 

STRT02 

SUBNOD 

SURFP_FSW 

T 

TC 

TDEPTH 

TDSAT 

TDTMAX 

TEXP 

TIME 

TLIN 

TFS 

TIME 

TMIN 

TRG 

TNODE 

TSAT 

TSEC 

TTIS 

TT 

TTO 

TT_SUM 

TTSTC_SUM 

TXOVER 


UFEET 

ULBL 

UMETER 

UNITS 

ZT 


RATE  OF  INSPIRED  02  TENSION  CHANGE  IN  CURRENT 
PROFILE  SUB-SEGMENT  (FSW/MIN) 

END  SURFACE  OR  SHALLOW  INTERVAL  TIME  PER  END- 

INTERVAL  REPETITIVE  GROUP  ARRAY 

TIME  (MIN)  REMAINING  IN  02  OR  AIR  BREATHING 

CYCLE  DURING  02 /AIR  CYCLING  WHEN  OXYD= . TRUE . 

SATURATION-DESATURATION  HALFTIME  RATIO  ARRAY 

BOOLEAN.  TRUE:  START  02 /AIR  CYCLING 

BOOLEAN  SET  WITH  "SD"  OPTION  IN  PCF  FILE. 

ABSOLUTE  SURFACE  PRESSURE  IN  FSW  UNITS 

ELAPSED  TIME  DIVE  PROFILE  ARRAY 

TIME  CHANGE  BETWEEN  PROFILE  SUB-SEGMENT  END 

NODES  (MIN) 

TARGET  ASCENT  DEPTH  (FSW  OR  MSW) 

STOP  TIME  FOR  DESATURATING  TISSUE  (MIN) 

MAXIMUM  PERMISSABLE  TOTAL  DIVE  TIME  (MIN) 

TIME  IN  A  PROFILE  SUB-SEGMENT  DURING  WHICH  GAS 
EXCHANGE  IS  EXPONENTIAL  (MIN) 

STOP  TIME  (MIN) 

TIME  IN  A  PROFILE  SUB-SEGMENT  DURING  WHICH  GAS 
EXCHANGE  IS  LINEAR  (MIN) 

ASCENT  TIME  TO  FIRST  STOP  (MIN) 

TIME  AT  CURRENT  DEPTH  (MIN) 

WHOLE  MINUTES  PORTION  OF  TIMES 

TIME  TO  REACH  A  GIVEN  REPETITIVE  GROUP 

DESIGNATOR  (MIN) 

TIME  OF  MINIMUM  OR  MAXIMUM  TISSUE  INERT  GAS 
TENSION  IN  A  PROFILE  SUB-SEGMENT  (MIN) 

STOP  TIME  FOR  SATURATING  TISSUE  (MIN) 

WHOLE  SECONDS  PORTION  OF  TIME 

BOOLEAN.  TRUE:  INCLUDE  TRAVEL  TIME  IN  STOP  TIME 
ASCENT  TRAVEL  TIME  TO  BE  INCLUDED  IN  NEXT  STOP 
TIME  (MIN) 

PREVIOUS  TT  VALUE  (MIN) 

CUMULATIVE  ASCENT  TRAVEL  TIME  NOT  YET  INCLUDED 
IN  TOTAL  DECOMPRESSION  TIME  (MIN) 

CUMULATIVE  ASCENT  TRAVEL  TIME  INCLUDED  IN  STOP 
TIME  (MIN) 

TIME  IN  A  PROFILE  SUB-SEGMENT  AT  WHICH 
CROSSOVER  TO  OTHER  GAS  TRANSPORT  MODE  OCCURS 
(MIN) 

FEET  LABEL 

DEPTH  UNITS  LABEL  FOR  PRINTOUT 
METERS  LABEL 

DEPTH  UNITS  INDICATOR  (1=FSW;  2=MSW) 

DIVE  PROFILE  CUMULATIVE  ELAPSED  TIME  ARRAY 
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7.  Appendix  B.  Planner  Coordinate  Format  (PCF)  Description 


The  Planner  Coordinate  Format  (PCF)  was  developed  to  support 
coding  of  incomplete  profiles  with  instructions  for  completing 
them.b  The  format  structure  reflects  the  usual  diving  practice  of 
effecting  breathing  gas  switches,  when  they  occur,  on  arrival  at 
stops  during  either  compression  or  decompression.  Complete 
profiles  can  also  be  coded  in  this  format  with  invocation  of  the 
global  GDX  directive  and  use  otherwise  of  only  the  TX  (GTX) ,  TS 
(GTS),  TZ  (GTZ) ,  and  AB  directives.  Enhancements  to  the  PCF 
developed  to  support  completion  of  profiles  with  probabilistic 
models  of  decompression  sickness  incidence  and  time  of  occurrence 
are  not  included  in  this  description. 

In  the  following,  code  elements  in  square  brackets  are  optional: 
the  brackets  are  not  part  of  the  element. 

The  first  two  lines  of  each  PCF  file  give  global  information  about 
how  all  subsequent  profiles  in  the  file  are  to  be  processed.  The 
remaining  lines  of  the  PCF  file  give  information  about  one  or  more 
separate  profiles  and  how  each  is  to  be  completed. 

First  two  lines  in  file  (in  1-2) : 

model  parameterization  filename 

pressure  unit  index,  stop  depth  increment,  global  options 

where  global  options  include:  [, TABLE] [,NAB] [,NBL] [,NPB] 

[ ,  GTX]  [ ,  GTS  ]  [ ,  GTZ  ]  [,GDX]  [,GLS##]  [,GOD]  [,GIW##] 

[ , CNDSDR] 

First  two  lines  of  profile  (in  1-2) : 

Freeform  comments 

FI,  PI,  F2 ,  P2 ,  F3 ,  P3 ,  F4,  P4,  ... 

Second  two  lines  of  profile  (in  3-4) : 

Dl,  D2 ,  RATE 

IGAS,  OPTIONS 


b  An  early  version  of  this  format  appeared  in  E.  D.  Thalmann,  Computer  Algorithms  Used  in  Computing  the 
MK  15/16  Constant  0.7  ATA  Oxygen  Partial  Pressure  Decompression  Tables ,  NEDU  TR  1-83,  Panama  City, 
FL,  Jan  1983. 
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These  profile  lines  specify  the  initial  conditions  for  the 
profile  and  the  first  stage  in  the  profile.  By  default,  the 
profile  is  started  from  initial  saturation  on  air  at  the 
pressure  D1  specified  in  this  record,  with  reference  surface 
pressure  at  sea  level.  The  first  stage  in  the  profile  is 
then  to  travel  to  D2  on  IGAS  at  the  specified  RATE. 

Each  remaining  pair  of  lines  in  the  profile  —  (in  5-6) ,  (in  7- 
8),  etc.  —  contains  entries  in  the  following  format: 

TIME,  D2 ,  RATE 
IGAS,  OPTIONS 

The  entries  describe  the  PRESSURE,  respired  gas,  and  TIME 
properties  —  or  provide  information  for  calculation  of  these 
properties  —  between  three  nodes  in  a  dive  profile  as 
schematized  in  Figure  B.l: 


Pressure 


D1  (=D2  from  node  2  of  previous  record  or  initial 
conditions  if  present  record  is  the  first  in  the  profile) 

*  RATE  (as  specified) 


D2 


Time 


Node  # : 


H -  TIME  - * 

U -  IGAS  - 


I _  *  Node  0  is  node  2  of  previous  record  or 

as  defined  by  initial  conditions  if  the 
present  record  is  the  first  in  the  profile. 

Figure  B.l.  Schematization  of  nodes  (filled  circles)  generated  from  a  two-line  *.pcf 
record.  Slashed  lines  indicate  that  some  information  about  pressure,  respired  gas, 
and  time  in  these  internodal  regions  may  be  calculated  in  accord  with  instructions  in 
the  OPTIONS  list.  Such  calculation  may  result  in  insertion  of  one  or  more  additional 
nodes  between  the  three  illustrated  nodes. 
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Each  two-line  record  can  be  expressed  in  the  following  words: 
"After  spending  TIME  at  Dl,  travel  to  D2  at  RATE  —  all  on 
I  GAS.  " 

Note  that  the  PCF  record  contains  no  explicit  information 
about  the  node  2  time  (Figure  B.l) .  This  time  is  computed 
with  D2  and  the  specified  RATE  by  default  methods,  or  as 
specified  by  instructions  for  completion  of  the  Dl  ->  D2 
transition  in  the  OPTIONS  list.  Available  options  are 
described  in  a  subsequent  section  of  this  description. 

If  D2  <  Dl,  the  default  is  to  calculate  and  insert 
decompression  stops  between  Dl  and  D2  for  decompression  on 
IGAS  from  Dl  to  D2 .  After  insertion  of  intervening  stops, 
default  processing  of  a  D2  <  Dl  *.pcf  record  completes  with 
arrival  at  D2 .  If  an  LS  option  is  specified,  the  record 
completes  after  insertion  of  any  time  required  at  D2  before 
ascent  can  continue  to  D2  -  DINC.  Of  course,  if  a  no-stop 
decompression  from  Dl  to  D2  can  be  effected  without  violation 
of  maximum  permissible  tissue  tensions  at  any  intervening 
depth  or  D2,  decompression  from  Dl  to  D2  at  RATE  is  effected 
without  stops. 

The  last  pair  of  lines  in  the  profile  must  include  the  FN  option 
to  indicate  that  the  pair  constitutes  the  last  instruction  for 
the  profile. 

Successive  profiles  in  the  PCF  file  can,  but  need  not,  be 
separated  by  one  or  more  blank  lines. 


Field  Descriptions 


model  parameterization  filename  --  Filename  of  model 

parameterization  file  to  use  when  processing  the  PCF 
file.  The  model  parameterization  input  file  (MPIF) 
and  its  associated  DCA_INIFL  (DCA_params_*MPIF* . ini) 
are  each  sought  first  in  any  path  specified  with  the 
MPIF  filename  in  the  PCF  file,  then  in  the  PCF  file 
path,  and  finally  in  the  home  directory  of  the 
running  instance  of  the  program. 

pressure  unit  index  --  One-  or  two-character  integer  entry  as  in 
the  Augmented  NMRI  Standard  format. 
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stop  depth  increment  --  Floating  point  or  integer  numerical  entry 
giving  the  pressure  differential  between 
decompression  stops  in  units  prescribed  by  the 
preceding  default  pressure  unit  index. 

global  options 

CNDSDR  Optional  global  switch  for  conditional  use  of  SDR^l . 0 
values  in  MPTT  table.  [SDR(I)  is  the 

saturation/desaturation  time  constant  ratio  in  these 
applications.]  When  placed  in  the  global  options 
list,  CNDSDR  causes  SDR (I)  td . 0  values  from  the  MPTT 
table  to  be  used  only  when  inspired  F02  exceeds  the 
value  of  CNDSDR_F02  set  from  the  DCA_Params.ini  file. 
SDR (I) =1.0  is  used  for  all  I  compartments  otherwise. 
CNDSDR_F02  is  defaulted  to  0.0  if  not  found  in 
DCA_Params.ini,  or  if  the  CNDSDR  switch  does  not 
appear  in  the  PCF  globals  line,  causing  uniform  use 
of  the  SDR (I)  values  as  set  in  the  MPTT  table. 

GTX,  GTS,  GTZ ,  GDX,  GOD,  --  Optional  global  switches  to 

effect  TX,  TS,  TZ,  DX,  and/or  OD  OPTIONS  (See  OPTIONS 
below)  throughout  all  profiles  in  the  file,  whether 
or  not  the  respective  switches,  TX,  TS,  TZ,  DX,  or  OD 
are  specified  in  individual  records. 

NOTE:  The  GTX  option  causes  only  preceding  descent 
travel  time  per  se  to  be  subtracted  from  isobaric 
segment  times  in  a  given  descent  from  surface.  Time 
at  any  intervening  isobaric  stages  or  "holds"  is  not 
subtracted.  The  effect  of  the  GTX  option  on  a  series 
of  records  for  a  given  descent  can  thus  be  different 
from  that  for  explicitly  specified  TX  options  in 
those  records.  The  GTX  option  overrides  TX  option 
behavior . 

GIW[##]  Optional  global  switch  to  specify  that  all 

decompressions  use  in-water  02  decompression  with  02 
breathing  starting  at  the  first  stop  at  or  shallower 
than  02CEIL  specified  in  the  DCA_Params.ini  file.  If 
the  optional  ##  entry  accompanies  the  GIW  switch, 
where  ##  is  an  integer  or  floating  point  numerical 
entry  of  up  to  8  characters,  the  first  02  breathing 
stop  in  each  computed  decompression  is  set  at  the 
first  decompression  stop  at  or  shallower  than  depth 
##,  irrespective  of  02CEIL  specified  in  the 
DCA_Params.ini  file.  ##  is  interpreted  in  the 
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default  depth  units  for  the  file  and  should  be  an 
integer  multiple  of  the  Stop  Depth  Increment 
specified  in  the  file  header.  NOTE:  The  gas  setting 
specified  by  the  prevailing  P*  or  F*  option  is  used 
until  02CEIL  (or  depth  =  ##)  is  reached  during 
ascent.  Air-breathing  breaks  are  inserted  throughout 
the  course  of  the  computed  02  breathing  times  as 
specified  by  the  02TIME  and  AIRTIME  settings  in  the 
DCA_Params.ini  file. 

GLS##  Optional  global  switch  to  set  the  last  or  shallowest 
allowed  stop  at  ##  in  all  computed  decompressions  in 
the  file,  where  ##  is  an  integer  or  floating  point 
numerical  entry  of  up  to  8  characters.  ##  is 
interpreted  in  the  default  depth  units  for  the  file 
and  should  be  an  integer  multiple  of  the  Stop  Depth 
Increment  specified  in  the  file  header. 

NAB  Optional  global  switch  that  suppresses  automatic 

insertion  of  air  breathing  breaks  in  02  breathing 
periods  during  02  decompression.  The  NAB  switch 
overrides  the  02TIME  setting  in  the  DCA_Params.ini 
f  ile . 


NBL 


NPB 


Optional  global  formatting  directive  to  specify  "NO¬ 
BLANK-LINES"  in  the  decompression  table  generated 
under  the  TABLE  switch. 


Optional  global 
PAGE  BREAKS"  in 
under  the  TABLE 


formatting  directive  to  specify  "NO 
the  decompression  table  generated 
switch . 


TABLE  Optional  switch  indicating  that  the  file  contains 

profiles  that,  when  completed,  can  be  output  in  the 
form  of  a  decompression  table.  Enables  the  PCF  file 
processor  to  generate  a  table  file  in  format 
approximating  that  used  in  the  US  Navy  Dive  Manual  as 
part  of  processor  output.  For  proper  function,  the 
PCF  file  should  be  arranged  with  profiles  in  order  of 
increasing  bottom  times  in  groups  of  increasing  dive 
depth . 


Dl,  D2  Floating  point  or  integer  entry  specifying  node 

pressures  in  units  prescribed  by  the  default  pressure 
unit  index  for  the  file. 
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FI,  F2,  ...,  Fn  Floating  point  numerical  entry  for  the  nth  inert  gas 

fraction  (or  Gas  Code)  breathed  in  the  profile, 
if  <  1.0,  Fn  =  N2  fraction  in  nth  respired  N2-O2  mix, 
dry . 

if  >  1.0,  Fn  =  Augmented  NMRI  Standard  format  Gas 
Code  for  nth  gas  breathed  in  the  profile.  Note  that 
this  is  the  only  means  by  which  respired  gases 
other  than  N2-O2  mixes  can  be  coded. 

IGAS  Alpha  "F"  or  "P"  followed  by  an  integer  pointing  to  a 

gas  description  in  the  2nd  line  of  the  PCF  file. 

Entry  indicates  the  gas  mix  breathed  during  the 
stages  described  in  the  current  profile  record. 

E.g.,  "F3"  indicates  that  the  third  gas  fraction  (or 
Gas  Code)  given  in  line  2  is  breathed  during  the 
current  profile  stages. 

OPTIONS  Series  of  upper  case  alpha  and  numeric  entries 

indicating  action  to  take  in  processing  the  present 
profile  record.  Options  may,  but  need  not,  be 
separated  by  commas  or  spaces.  More  than  one  option 
may  be  specified,  as  appropriate,  in  any  order. 

PI,  P2,  ...,  Pn  Floating  point  numerical  entry  for  the  nth  O2 
partial  pressure  breathed  in  the  profile,  wet 
(after  adjustment  for  water  vapor  pressure) . 

RATE  Floating  point  or  integer  numerical  entry  specifying 

the  rate,  in  UNITS/min,  for  the  D1  to  D2  transition 
(compression  or  decompression)  specified  in  the 
record.  The  UNITS  used  to  interpret  the  entry  are 
those  used  to  interpret  the  D2  entry  in  the  record. 

If  D2  =  Dl,  the  RATE  entry  has  no  effect. 

TIME  Floating  point  or  integer  entry  specifying  duration 

of  residence  at  Dl  in  minutes. 


Options 

AB  Insert  air  breaks  into  residence  on  IGAS  (FI_2)  at  Dl 

according  to  O2/AIR  cycling  specified  in  the  DCA_Params.ini 
file.  This  option  is  exercised  only  on  stages  for  which  Dl 
>  Surface  Pressure.  If  AB  is  specified  in  a  record  for 
which  a  decompression  is  computed,  O2/AIR  cycling  is 
continued  through  the  decompression  with  an  AIRTIME  minute 
air  break  after  each  02TIME  minutes  on  IGAS,  unless  an  OD 
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option  accompanies  the  AB  option.  In  the  latter  cases,  the 
OD  option  causes  gas  with  02TIME_F02  specified  in  the 
DCA_Params.ini  file  to  be  breathed  during  each  02TIME  period 
during  the  decompression,  but  not  during  preceding  residence 
at  D1 . 

DX  Decompress  from  D1  ->  D2  without  staging.  The  DX  option 

overrides  default  decompression  optimization  and  effects  no¬ 
stop  decompression  from  D1  to  D2 . 

FN  (or  FL)  Last  PCF  record  for  profile. 

IW  Invokes  the  OD  option  with  specific  settings  for  in-water  02 
decompression;  02  at  02TIME_F02  starting  at  02CEIL  specified 
in  DCA_Params.ini  file.  If  the  NAB  global  option  has  not 
been  set  at  the  beginning  of  the  file,  air-breathing  breaks 
are  inserted  throughout  the  course  of  the  computed  02 
breathing  time  of  the  decompression  as  specified  by  the 
02TIME  and  AIRTIME  settings  in  the  DCA_Params.ini  file. 

NOTE:  The  gas  setting  specified  by  the  leading  P*  or  F* 
option  in  the  record  is  used  until  02CEIL  is  reached  during 
the  ascent. 

LS  Compute  decompression  to  D2  with  all  stops  allowed  as  by 

default,  or  as  specified  by  a  GLS##  directive,  and  include 
any  required  time  at  D2  in  the  final  profile.  A  local  LS 
option  overrides  any  setting  from  a  global  GLS##  setting  for 
the  decompression  specified.  The  LS  option  has  no  effect  in 
the  presence  of  the  DX  (or  GDX)  option,  if  D2  >=  Dl,  or  if 
D2  =  Surface  Pressure. 

OD  Compute  decompression  using  oxygen  breathing.  If  the  NAB 

global  option  has  not  been  set  at  the  beginning  of  the  file, 
air-breathing  breaks  are  inserted  throughout  the  course  of 
the  computed  decompression  as  specified  by  the  02TIME  and 
AIRTIME  settings  in  the  DCA_Params.ini  file. 

SD  Generate  a  row  in  a  Sub-NoD  table  consisting  of  the  no-stop 
limit  and  the  no-stop  times  required  to  reach  the  top  of 
each  surfacing  repetitive  group  for  depth  =  D2 . 

TS  Include  travel  time  to  each  stop  depth  in  computed  stop 

times  if  the  stop  is  not  the  first  resulting  from  execution 
of  the  instruction.  This  option  has  no  effect  with  a  DX  or 
GDX  option,  or  if  D2  >=  Dl . 
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TX  Interpret  TIME  as  Bottom  Time  (subtract  preceding  travel 
time  from  coded  TIME  value) .  This  option  only  affects 
isobaric  stages  for  which  D2  >  D1 .  All  preceding  time  since 
leaving  surface  or  since  start  of  the  last  isobaric  segment 
adjusted  with  a  TX  option  after  leaving  surface,  including 
any  "hold"  times  at  isobaric  segments,  is  subtracted  from 
the  time  specified  in  the  record.  A  prevailing  global  GTX 
option  overrides  any  TX  option  specifications. 

TZ  Reduce  specified  surface  interval  by  total  decompression 

time  of  preceding  dive.  The  TZ  option  is  implemented  after 
completion  of  any  specified  gas  switch  on  arrival  at  surface 
pressure,  and  only  on  the  first  segment  explicitly  coded  for 
a  given  surface  interval.  This  OPTION  has  no  effect  if  D2  ^ 
D1  ^  Surface  Pressure. 
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